Vortex Lattice Structural Transitions: a Ginzburg-Landau Model Approach 
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We analyze the rhombic to square vortex lattice phase transition in anisotropic superconductors 
using a variant of Ginzburg-Landau (GL) theory. The mean-field phase diagram is determined to 
second order in the anisotropy parameter, and shows a reorientation transition of the square vortex 
lattice with respect to the crystal lattice. We then derive the long-wavelength elastic moduli of the 
lattices, and use them to show that thermal fluctuations produce a reentrant rhombic to square 
lattice transition line, similar to recent studies which used a nonlocal London model. 
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PACS numbers: 74.25.Qt,74.40.-|-k,74.70.Dd 

Vortex lattices in superconductors with four-fold aniso- 
tropy can exhibit a rhombic to square structural phase 
transition due to a competition between the anisotropy 
(favoring a square lattice) and the repulsive vortex- vortex 
interaction (favoring a triangular lattice) . This transition 
has been observed in the borocarbide family of supercon- 
ductors i?Ni2B2C (with R a rare earth element), whose 
anisotropy is due to lattice effects Q, 0, S S Pli and 
more recently in the overdoped cuprate superconductor 
Lai.83Sro.i7Cu04+5 whose anisotropy is due to un- 
conventional superconducting order. Recent small angle 
neutron scattering (SANS) data for LuNi2B2C JJ show 
the transition line in the H — T plane curving upward 
and avoiding Hc2, contrasting with mean- field predic- 
tions and underscoring the importance of thermal fluc- 
tuations even in this low temperature superconductor. 
Some features of the data can be understood using a 
nonlocal London model developed by Kogan and collab- 
orators 0il3; however, the London model is ostensibly 
only applicable at low magnetic inductions, leaving as an 
open question the behavior close to Hc2- An alternate ex- 
planation for the reentrant behavior of the vortex lattice 
transition is also available in the literature . 

In this Letter we study the effect of thermal fluctua- 
tions on the rhombic to square vortex lattice transition 
using an anisotropic GL model. We first determine the 
mean-field phase diagram to second order in the aniso- 
tropy parameter, and find that even at this level GL the- 
ory predicts a novel reorientation transition of the square 
vortex lattice with respect to the crystal lattice, an effect 
overlooked in previous GL studies. We then derive the 
(nondispersive) elastic moduli of the lattices using meth- 
ods first pioneered by Brandt for isotropic supercon- 
ductors. Finally, these moduli are used to determine the 
fluctuation renormalized phase boundary in the H — T 
plane. For parameters typical of the borocarbides we find 
a reentrant phase boundary, in agreement with the non- 
local London model results of Ref. ■ As GL theory is 
expected to be valid near Hc2, our results nicely comple- 
ment and extend the nonlocal London model results. 

Accounting for anisotropy, the linearized first GL equa- 



tion (in dimensionless units) is 

tfV - V' - e(nl - n2)2v> = 0, 



(1) 



where 11 = [W /in — A) is the covariant derivative and 
e is the dimensionless anisotropy parameter. For the 
cuprates the anisotropic term can be derived rigorously 
from the two field s and d-wave GL free energy [ij, [iJl 
by integrating out the s-wave term [l^ . For the conven- 
tional, but highly anisotropic, borocarbides it is a phe- 
nomenological term that will be sufficient to capture all 
the interesting physics. The symmetry breaking proper- 
ties of the anisotropic term dictate the introduction of 
the orientation angle (f) of the vortex lattice relative to 
the underlying crystal lattice — this term provides an ef- 
fective coupling of the vortex and crystal lattices, and 
makes structural transitions possible. 

We have generalized the virial theorem of supercon- 
ductivity jlSl Il6j to include the anisotropic term in the 
extended GL theory 0|, with the result 
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where the brackets () denote spatial averaging. 

Using the generalizations of the Abrikosov identities 
derived with the help of Eq. |(5J , the free energy can be 
cast into the compact form 
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where (Sa is the usual Abrikosov parameter and F is the 
correction to the isotropic result, 
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Re{r{Til-I^lf^|;)^ 



+ (2^2-1) 



Throughout this paper we will use the convenient param- 
eterization p + ia = (b/a) exp(z6'), where b and a are the 



2 



magnitudes of the two basis vectors of the vortex lattice 
and 9 is the apex angle. 

Equation ^ was solved perturbatively to second or- 
der in the anisotropy par ameter e, extending previous 
first order results jll. Il9l Eoj . The second order cal- 
culation is necessary to justify the use of perturbation 
theory and at the same time to quantify the accuracy of 
the first order correction. More importantly, we uncover 
new and unexpected physics at second order: a reorien- 
tation transition of the square vortex lattice with respect 
to the crystal lattice. The quantities Pa and T are given 
in terms of rapidly converging sums involving p, cr, and 



r = 8e62|/3A(2K2 - 1)[H- 3(efe)] + /3o 

+ /3i[l + 6(e6)] + 2(e6)[7i+72] 
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where the functions /?o, 7i, 72 are periodic in p with 
period 1, and are given by 
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where A„™ = ^(P-i^P(n'-ra^)^-2^<y(n'^-m^) ^ rpj^g pj.jj^g 

denotes a double sum: one over n, m and a second with n 
and m replaced by (n -I- 1/2) and (m -f- 1/2) respectively. 
We perform a numerical minimization of the free energy 
(PJ with respect to p, a, and <^ to obtain the mean-field 
vortex lattice structure as a function of the field and tem- 
perature. We find a continuous transition from a rhom- 
bic to square phase as e is increased. The second order 
correction only shifts the transition point by about 1% 
relative to the first order calculation, contrary to previ- 
ous claims [2^ which neglected the lifting of the orienta- 
tion degeneracy by the anisotropic term. Our numerical 
investigations show unambiguously that a consistent ex- 
ploration of the relevant phase space and an accurate 



determination of the global minima of the GL free en- 
ergy lO demands the inclusion of the the orientational 
degree of freedom 0. After the structural phase tran- 
sition occurs the square vortex lattice is oriented along 
the [110] direction, and remains in this preferred orien- 
tation until a field hr is reached. At that point another 
continuous phase transition takes place, with the vortex 
lattice changing its orientation while retaining its square 
shape. At a given b > br, two equivalent configurations 
exist: one with orientation angle (j) > 7r/4, and one with 
<f>' = n /2 — (j). This reorientation transition occurs only 
at second order in the anisotropy e. 

Having determined the mean-field phase diagram of 
the system, we now turn to the elastic properties of the 
lattices. Assuming that the displacement field u is small 
compared to typical lattice spacings, we calculate the 
elastic moduli from the GL free energy 10) in the in- 
compressible limit, V • u = 0. For simplicity wc have 
focused on straight and parallel vortices; we do not ex- 
pect the structural transition to significantly affect the 
tilt modulus C44. The relevant moduli are essentially non- 
dispersive allowing us to consider uniform deforma- 
tions. Now there are many ways to choose the elastic 
moduli from the independent components of the elastic 
tensor, and we will follow the lead of Miranovic and Ko- 
gan |22i | and use as the four independent elastic moduli 
for a rhombic lattice the quantities Csq, cesx, cgey, and c^, 
corresponding to four particular deformations. The so- 
called "squash" deformation is responsible for the trans- 
formation of the rhombic to a square vortex lattice, and 
the corresponding squash modulus Csq vanishes at the 
transition point. The other relevant deformations are 
the two pure shears in x and y directions (which deter- 
mine cesx, ceey) and the rotation of the lattice about the 
direction of the magnetic field (which determines c, ). 

All elastic moduli are obtained from the second deriva- 
tive of the GL free energy © . The elastic energy can be 
written as 



(11) 



where c(/i) is the elastic modulus for each particular de- 
formation. Each elastic modulus is obtained from 



+ 2Cpap^ap + 2Cp^pp,(f)p_ + 2Ccr(pap4>p, 



(12) 



where Cy are the terms originating from the derivatives 
of the GL free energy (j2Jl 



= (2~<5,,) <i>(6^-F)/3A,,,+e 



+ (3a,jT^, + r/3A,y - 2^TPa,^Pa,j) 



(13) 



The derivatives in the right hand side of Eq. I|13|l are with 
respect to i,j = p, a, cj), ^ = {2k^ - 1)/[{2k'^ ~ 1)/3a + 1], 
and e = (1 - 6)V[(2k2 _ l)/5^ + l]2. 
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The calculated elastic moduli for e = 0.11 and k = 10 
(appropriate for the borocarbides) are shown in Fig. ^ 
The transition from a rhombic to a square vortex lattice 
is signaled by the vanishing of the squash modulus. At 
the same point the two shear moduli merge into one be- 
cause of the higher symmetry of the square phase. At a 
still higher field we observe the reorientation transition, 
which occurs when = 0. This instability of the vortex 
lattice with respect to rotations also exists in the nonlocal 
London model "2^. There are experimental indications 
for this instability, although its investigation has been 
hampered so far by the high fields required for its obser- 
vation [2^ . This reorientation transition is sensitive to 
the value of k, and from our results we predict that this 
transition will not occur in the cuprates (k ^ 100), as 
the transition point where Cr = moves above Hc2- 

We also note that anisotropic elasticity with a softer 
shear modulus for shearing along the sides of the square 
lattice than along the diagonals emerges naturally in this 
model. This response of the vortex lattice is measured by 
the ratio C66a:/c66,f 1 which has a maximum at the tran- 
sition as can be seen in Fig. |21 It has been suggested 
[23! that this behavior explains the anisotropic orienta- 
tional long range order observed in decoration patterns in 
LuNi2B2C, which manifests itself as a significant differ- 
ence in the correlation lengths along the [110] and [110] 
directions. 

We next turn to the effects of thermal fluctuations on 
the vortex lattice structure. Following Gurevich and Ko- 
gan [lll |. we will make the simplifying assumption that 
harmonic fluctuations renormalize the three variables p, 
fj, and (f), which characterize the structure of the vortex 
lattice, with the result that 

p = p[l-a(2-p)^], (14) 
CT==cr[l-cr(2-p)^], (15) 

(j) = (j) + a cot (j) u'^ , (16) 

where overbars denote thermal averages. The rhombic 
to square transition is then determined as the locus of 
points in the H-T plane where c^q — 0. For that pur- 
pose, the "squash" elastic modulus Csq is recalculated 
from Eq. H12|l using the renormalized variables p, W, (p. 
This is not performed self-consistently, as we do not take 
into account the effects of the fluctuations on the elas- 
tic moduli themselves. This approach is justified due 
to the fact that the mean squared displacement of the 
vortex v? is finite at Hd{T) 01 ■ There are two com- 
peting parameters, the strength of the thermal fluctu- 
ations X the anisotropy parameter e which is re- 
lated to the nonlocality parameter pcK of the nonlocal 
London model as e = {pcK/^af (1 — t^)/12. To com- 
pare with the nonlocal London model results, we take 
Inl A(r) = Ao/yr^ and jjT) = ^o/TT^, where 
t = T/Tc, so that X = XoVvl^^- The dimensionless 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



b 

FIG. 1: The elastic moduli versus the reduced field h for 
e = 0.11. In the inset we show the region where the reori- 
entation transition happens, as is signaled by the vanishing 
of the rotation modulus Cr- 
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FIG. 2: The ratio ceei/cee.-j versus the reduced field 6. 

mean squared displacement is then 

^ ^ Xo^ ht 

y(r^[(l-6)3 1n(l-H^)]^' 

with ?7 ~ 3 for LuNi2B2C 11]. The results of the numer- 
ical solution are shown in Fig. |31 for the same param- 
eter values used by Gurevich and Kogan namely, 
PCK — 2.5^0 and xo = 0.0064. With fluctuations, the 
phase boundary curves backward and displays the "reen- 
trant" behavior that was also obtained by Gurevich and 
Kogan. The physical explanation of this feature is sim- 
ple: thermal fluctuations tend to wash out effects of the 
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Transition lines for p=2.5t , X =0.0064 



— GL with fluctuations 

— GL without fluctuations 

— Non-local London model 




T/T 



FIG. 3: The phase boundaries in the H-T plane separating the 
square and rhombic vortex lattice obtained from GL theory 
and the nonlocal London model. The line intercepting Hc2{t) 
is the mean field GL result. 



anisotropy, so that the preferred symmetry of the lat- 
tice is rhombic. Note that the abrupt termination of the 
transition line for the GL theory is an artifact of the per- 
turbative treatment of the anisotropy; at that point the 
correction term F ceases to be small and the extended GL 
theory becomes unstable. While this shortcoming of the 
model could be remedied by including higher order terms 
in the GL free energy, the overall qualitative features of 
the phase diagram are expected to remain unchanged. 
One such choice would be an additional isotropic term. 
Nevertheless, the complications introduced in the anal- 
ysis by the unnecessary extra parameter, would clearly 
outweigh any potential benefits. 

We have used the GL theory with an anisotropic, sym- 
metry breaking term to study the evolution of the vortex 
lattice and particularly the structural phase transition 
and its connection with the elastic response of the vor- 
tex lattice. We have shown that to obtain the correct 



phase diagram it is necessary to calculate properties up 
to second order in the anisotropy; at this order a rota- 
tional instability of the square vortex lattice appears. We 
predict that this instability should be observable in the 
borocarbides but not in the cuprate superconductors. In- 
cluding thermal fluctuations in the model we were able to 
show that the rhombic to square transition line exhibits 
reentrant behavior, similar to what is obtained within 
the nonlocal London model, and in qualitative agreement 
with SANS studies of LuNiaBaC. 

We would like to thank V. G. Kogan, A. Gurevich and 
P. J. Hirschfeld for useful discussions. This work was 
supported by NSF DMR-9978547. 



[1 

[2; 
[3; 
[4: 

[5 
[G 
[7] 
[8] 

[9: 

[lo; 

[11 

[12 

[13 
[14 

[15 

[16 
[17 
[18 
[19 
[20 

[21 
[22; 

[23] 



M. R. Eskildsen et al, Phys. Rev. Lett. 86, 5148 (2001). 
M. R. Eskildsen et al., Phys. Rev. Lett. 78, 1968 (1997). 
Y. De Wilde et al, Phys. Rev. Lett. 78, 4273 (1997). 
Sakata et al, Phys. Rev. Lett. 84, 1583 (2000). 
Ya. Vinnikov et al, Phys. Rev. B 64, 024504 (2001). 
Gilardi et al, Phys. Rev. Lett. 88, 217003 (2002). 
G. Kogan et al, Phys. Rev. B 54, 12386 (1996). 

G. Kogan et al, Phys. Rev. B 55, R8693 (1997). 
Nakai et al, Phys. Rev. Lett. 89, 237004 (2002). 

H. Brandt, Rep. Prog. Phys. 58, 1465 (1995). 
Gurevich and V. G. Kogan, Phys. Rev. Lett. 87, 

177009 (2001). 

M. Franz, I. Affleck, and M. H. S. Amin, Phys. Rev. Lett. 
79, 1555 (1997). 

R. Heeb et al, Phys. Rev. B 54, 9385 (1996). 

I. Affleck, M. Franz, and M. H. S. Amin, Phys. Rev. B 

55, R704 (1997). 

M. M. Doria, J. E. Gubernatis, and D. Rainer, Phys. 
Rev. B 39, 9573 (1989). 

U. Klein and B. Pottinger, Phys. Rev. B 44, 7704 (1991). 

A. D. Klironomos and A. T. Dorsey (unpublished). 
D. Chang et al, Phys. Rev. B 57, 7955 (1997). 

D. Chang et al, Phys. Rev. Lett. 80, 145 (1998). 

B. Rosenstein and A. Knigavko, Phys. Rev. Lett. 83, 844 
(1999). 

K. Park and D. A. Huse, Phys. Rev. B 58, 9427 (1998). 
P. Miranovic and V. G. Kogan, Phys. Rev. Lett. 87, 
137002 (2001). 

L. Ya. Vinnikov et al, Phys. Rev. B 64, 220508 (2001). 



